Spectral resolution in hyperbolic orbifolds, 
quantum chaos, and cosmology 



H. Then 



Abstract 



We present a few subjects from physics that have one in common: the 
spectral resolution of the Laplacian. 



1 Introduction 

If we look to nature, we observe dynamics and structure formation in various as- 
pects. Several theories exist that explain the growth of structure quantitatively. 
On the smallest scale, it is quantum mechanics that governs the dynamics, 
whereas on the largest scale the evolution of our universe follows the Einstein 
field equations. 

For all scales and almost any kind of observation we have equations at hand 
that quantify our findings. Many of these equations contain the spectral reso- 
lution of the Laplacian. Just think of the Schrddinger equation. In absence of 
forces or when the potential can be transformed into the metric, the station- 
ary Schrddinger equation reduces to the eigenvalue equation of the Laplacian 
subject to some boundary conditions. 

Spectral resolution has important consequences on many topics. Focusing 
on a few particular subjects, we start with an example from thermodynamics 
and demonstrate the importance of the spectral density. According to Weyl's 
law, the spectral density in leading order does not depend on the shape of 
the boundaries. In consequence, the properties of thermodynamic systems are 
universal. 

Less universal is the distribution of the eigenvalues. A central subject of 
quantum chaos is to classify the distribution of eigenvalues in dependence of 
whether the corresponding classical system is chaotic. In an explicit example, 
we compute the spectrum of the Laplacian numerically and confirm a conjecture 
of arithmetic quantum chaos. 

Finally, we use the eigenvalues and eigenfunctions of the Laplacian to com- 
pute the temperature fluctuations in the cosmic microwave background. 
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2 Thermodynamics 

The thermodynamic properties of an ideal gas can be deduced from the loga- 
rithm of the partition function, 

InZ = -r/^ln(l - 7 ? zc-' 3e "), 

a 

where a stands for the quantum numbers of a single particle with energy e a , 
z = e^ is the fugacity, j3 = 7™ is the inverse temperature, and 77 is the statistic 
parameter that distinguishes between bosons, 77 = 1, and fermions, r\ = — 1. 
Introducing the spectral density 

d(e) := y^S(e - e a ) 

a. 

the partition function can be expressed by an integral 

/oo 
d(e) ln(l - rjze- ffe ) de 
-00 

with the advantage that it is easier to compute an integral than a sum analyti- 
cally. 

Example 1 (The ideal gas inside a cube) The quantum microstate of an 
ideal gas in thermal equilibrium is subject to the stationary Schrddinger equation 



with Dirichlet boundary conditions at x,y,z = and X,y,Z = L. Since the 
Hilbert space of an ideal gas separates, the Hamiltonian H can be expressed by 
the sum of single-particle Hamiltonians hj = — w—Aj, j = 1, . . . , N, where the 
Laplacian Aj acts only on the coordinates of the j-th particle. Solving hip = sip 
yields the single-particle wave functions 

1 n ■ (Px x \ n ■ (Pvv\ n ■ (p* z 

and the single- particle eigenvalues 



Pl+P 2 y+P 2 z _ h 2 n 2 2 2 2 

lz — n ~ ^ TO \ n x ' n V "T 11 z) 



2m 2mL 2X x v 

with the quantum numbers n x ,n y ,n z S IN. 

In the thermodynamic limit, V — ■> oo, the spectrum becomes continuous. 
Consequently, the sum over the quantum numbers, 



E = EEE E 



n x n v n z m— — s 
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can be replaced by integrals 
resulting in 

j3„ 



— ^(2s + l) / d : 

4^ 2^ 3 y E 3 



(2tt/Y) 

where V — I? is the volume and s is the spin. 
Evaluating the spectral density 



P, 

R 3 



with 



yields 



d( £ )=^<5( £ - £Q )^^^(2.s + l) J d 3 p5(e-e a ) 



8(e - Ea )=5(e-£) = ^ (S(\p\ - sFhnS) + S(\p\ + V^e)) . 



(2m)i 

die) ~ (2s + 1) \ ' Vyfs for e > and die) = for e < 0. 
Inserting it into the partition function gives 

T/ o r°° t 

\nZ = -n{2s+\)-g-j=] dxv^Ml-^O, A: =V /^/3. 
Taylor expanding the integrand 

ln(l - V ze- X ) =-Y ^ ZC ~ X)n = V ze- X + -z 2 c~ 2x + 0(z 3 ), 
^-^ n 2 

n=l 

and evaluating the conditioning equation for the chemical potential, we have 

(2.+ !,£,(! + ^ + <V,). 
Hence, the thermal equation of state 

J:=~\nZ=-pV, 

P 

results in the virial expansion of the ideal gas 

77 TV A 3 , 



pV = NkTll 



2s+lVAV2 
Interested in the inner energie, we obtain 

dp 2 /3 2 F 



For further details and insight into thermodynamics, we refer the reader to any 
standard textbook on statistical physics, e.g. [15]. 

But what, if we ask for the ideal gas being inside a 3-sphere? Do we need 
to compute the entire example again, beginning with the spectral resolution of 
the Laplacian in a sphere? 

Theorem 1 (Weyl's law [29, 4]) If a quantum system is restricted to a finite 
volume V in D dimensions, its level counting function 

N(E) := #{ Ei \Ei<E}= f d(E) dE 

J — oo 

is asymptotically equal to 

N(E) ~ (2s + l)C D VE% 
in the semiclassical limit, E — > oo ; where 

r _ _i ( m \f 

D r(f + i)W j 

is a universal constant that only depends on the dimension D, but not on the 
specific shape of the boundary. 

Weyl's law is nice. It tells us that the results of thermodynamics are univer- 
sal, i.e. they are independent of the boundary shape. 



3 Quantum chaos 

In terms of classical mechanics a system is described by specifying the values of 
its coordinates and velocities, x%, ■ ■ ■ , Xf and x±, . . . ,if, where / is the number 
of degrees of freedom. 

If there exist / linear independent constants of motion that are in involution 
with each other, there is a set of canonical coordinates on the phase space, the 
action-angle variables. The action variables are constants of motion and the 
angle variables are the natural linear, periodic coordinates on the torus. The 
motion on the torus is linear in the angle variables, and the system is called 
classically integrable. 

If there do not exist / linear independent constants of motion, the system is 
called classical non-integrable or chaotic. The dynamics is non-linear and shows 
an exponential sensitivity to initial conditions. 

Being a limiting case of quantum mechanics, one might expect to see the 
properties of classical mechanics in quantum theory. But it is not this simple. 
The Schrodinger equation is linear and there is no exponential sensitivity to 
initial conditions. Quantum chaos relics on the behaviour of the correspond- 
ing classical system. A quantum system is called chaotic if and only if the 
corresponding classical system is non-integrable. 
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Central questions of quantum chaos concern the eigenvalue statistics and the 
distribution of eigenvalues in the semiclassical limit. 

We use the following assumptions: The quantum mechanical system is desym- 
metrised with respect to all its unitary symmetries, and whenever we examine 
the distribution of the eigenvalues, we regard them on the scale of the mean 
level spacings. Moreover, it is believed that after desymmetrisation a generic 
quantum Hamiltonian possesses no degenerate eigenvalues. 

Conjecture 1 (Berry, Tabor [6]) If the corresponding classical system is in- 
tegrable, the eigenvalues behave like independent random variables and the dis- 
tribution of the nearest-neighbour spacings is close to the Poisson distribution, 
i.e. there is no level repulsion. 

Conjecture 2 (Bohigas, Giannoni, Schmit [8, 9]) If the corresponding clas- 
sical system is chaotic, the eigenvalues are distributed like the eigenvalues of 
hermitian random matrices. The corresponding ensembles depend only on the 
symmetries of the system: 

• For chaotic systems without time-reversal invariance the distribution of 
the eigenvalues should be close to the distribution of the Gaussian Unitary 
Ensemble (GUE) which is characterised by a quadratic level repulsion. 

• For chaotic systems with time-reversal invariance and integer spin the 
distribution of the eigenvalues should be close to the distribution of the 
Gaussian Orthogonal Ensemble (GOE) which is characterised by a linear 
level repulsion. 

• For chaotic systems with time-reversal invariance and half-integer spin 
the distribution of the eigenvalues should be close to the distribution of the 
Gaussian Symplectic Ensemble ( GSE) which is characterised by a guartic 
level repulsion. 

These conjectures are very well confirmed by numerical calculations, but several 
exceptions are known. Here are two examples: 

Exception 1 The harmonic oscillator is classically integrable, but its spectrum 
is eguidistant. 

Exception 2 The geodesic motion on surfaces with constant negative curva- 
ture provides a prime example for classical chaos. In some cases, however, 
the nearest-neighbour distribution of the eigenvalues of the Laplacian on these 
surfaces appears to be Poissonian. 

Conjecture 3 (Arithmetic Quantum Chaos [7, 10]) On surfaces of con- 
stant negative curvature that are generated by arithmetic fundamental groups, 
the distribution of the eigenvalues of the guantum Hamiltonian are close to the 
Poisson distribution. Due to level clustering small spacings occur comparably 
often. 
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In the next sections, we compute numerically the eigenvalues and eigenfunc- 
tions of the Laplacian that describe the quantum mechanics of a point particle 
moving freely in the non-integrable three-dimensional hyperbolic space of con- 
stant negative curvature generated by the Picard group. The Picard group is 
arithmetic and we find that our results are in accordance with the conjecture of 
arithmetic quantum chaos. 

4 The modular surface 

For simplicity we first introduce the topology and geometry of the two-dimensional 
surface of constant negative curvature that is generated by the modular group 
[26] . It will then be easy to carry over to the three-dimensional space of constant 
negative curvature that is generated by the Picard group. 
The construction begins with the upper half-plane, 

H = {(x,y)eR 2 , y>0}, 

equipped with the hyperbolic metric of constant negative curvature 

, 2 dx 2 + dy 2 
ds = - . 

V 

A free particle on the upper half-plane moves along geodesies, which are 
straight lines and semicircles perpendicular to the x-axis, respectively. Express- 
ing a point (x, y) € TL as a complex number z = x + iy, all isometries of the 
hyperbolic metric are given by the group of linear fractional transformations, 

az + b 

zt-^^z^ -, a, o, c, a £ R, ad — bc=l, 

cz + d 

which is isomorphic to the group of matrices 

up to a common sign of the matrix entries, 

SL(2,R)/{±1} = PSL(2,R). 

In analogy to the concept of a fundamental cell in a regular lattice of a crystal 
we can introduce a fundamental domain of a discrete group T C PSL(2, R). 

Definition 1 A fundamental domain of the discrete group T is an open subset 
T C TL with the following conditions: The closure of T meets each orbit Tz = 
{72, 7 € r} at least once, T meets each orbit Tz at most once, and the boundary 
of T has Lebesgue measure zero. 
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If we choose the group T to be the modular group, 



T = PSL(2,Z), 
which is generated by a translation and an inversion, 

; j)=,», + i a» d (; 

the fundamental domain of standard shape is 

F={z = x + iy G H, -\< x <\> \ z \> 1 }- 

The isometric copies of the fundamental domain 7 J 7 , 7 G T, tesscllate the 
upper half-plane completely without any overlap or gap. Identifying the fun- 
damental domain T and parts of its boundary with all its isometric copies 
7JF, V7 G r, defines the topology to be the quotient space T\H. The quotient 
space T\H can also be thought of as the fundamental domain T with its faces 
glued according to the elements of the group T. 

Any function being defined on the upper half-plane that is invariant under 
linear fractional transformations, 

f(z) = f( 7 z) v 7 er, 

can be identified with a function living on the quotient space T\Tt and vice 
versa. 

With the hyperbolic metric the quotient space T\H inherits the structure of 
an orbifold. An orbifold locally looks like a manifold, with the exception that 
it is allowed to have elliptic fixed-points. 

The orbifold of the modular group has one parabolic and two elliptic fixed- 
points, 

a 1 -V3 

z = 100, z = 1, and z = — h 1 — . 

2 2 

The parabolic one fixes a cusp at z = ioo which is invariant under the parabolic 
element 

1 1 
1 

Hence, the orbifold of the modular group is non-compact. The volume element 
corresponding to the hyperbolic metric reads 

dxdy 
= — 5-, 

y 

such that the volume of the orbifold T\H is finite, 

voi(r\w) = |. 
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Scaling the units such that h = 1 and 2m = 1, the stationary Schrodinger 
equation which describes the quantum mechanics of a point particle moving 
freely in the orbifold T\H becomes 

(A + A)/(z) = 0, 

where the hyperbolic Laplacian is given by 



d 2 d 2 



k dx 2 dy 2 

and A is the scaled energy. We can relate the eigenvalue problem defined on the 
orbifold T\H to the eigenvalue problem defined on the upper-half space, with 
the eigenfunctions being subject to the automorphy condition relative to the 
discrete group T, 

f(r/z) = f(z) v 7 er. 

In order to avoid solutions that grow exponentially in the cusp, we impose the 
boundary condition 

f(z) = 0(y K ) for z^ioc, 

where k is some positive constant. 

The solutions of this eigenvalue problem can be identified with Maass wave- 
forms [17]. Numerically, they are most efficiently computed using Hejhal's al- 
gorithm [14]. 



5 The Picard surface 

In the three-dimensional case one considers the upper-half space, 

H = {(x ,x 1 ,y)eR 3 , y>0} 

equipped with the hyperbolic metric 

, n dx 2 , + dx 2 + dy 2 
ds 2 = — + — . 

y 

The geodesies of a particle moving freely in the upper half-space are straight 
lines and semicircles perpendicular to the xo-xi-planc, respectively. 

Expressing any point (xq, xi, y) G TL as a quaternion, z — xo + ixi + jy, with 
the multiplication defined by i 2 = —1, j 2 = —1, ij +ji = 0, all motions in the 
upper half-space are given by linear fractional transformations 

z i— > -fz = (az + b)(cz + d) , a, b,c,d G C, ad — be = 1. 

The group of these transformations is isomorphic to the group of matrices 

7 =(« ^)gSL(2,C) 
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up to a common sign of the matrix entries, 

SL(2,C)/{±1} = PSL(2,C). 

The motions provided by the elements of PSL(2, C) exhaust all orientation pre- 
serving isometries of the hyperbolic metric on 7i. 

We now choose the discrete group T C PSL(2, C) that is generated by the 
cosets of the following elements, 

1 l\ A i\ /0 -1 

o ij' v° V' V 1 

which yield two translations and one inversion, 

ZH>2+1, Z H z + i, Z t—¥ — z^ 1 . 

This group T is called the Picard group. The three motions generating T, 
together with the coset of the element 

i 

-i 

that is isomorphic to the symmetry 

z = x + jy i-> izi = -x + ]y, 
can be used to construct the fundamental domain of standard shape 

T = {z = x + ixi + jy G H, -- < x < -, < x x < -, \z\>l}. 

Identifying the faces of the fundamental domain according to the elements 
of the group T leads to a realisation of the quotient space T\H. 

With the hyperbolic metric the quotient space T\H inherits the structure of 
an orbifold that has one parabolic and four elliptic fixed-points, 

1 . [3 1 .1 . [T .1 



z = ioo, 2 = 1, z= — Hi - , z= — hi — Hi - , z = i — Hi - . 

2 1 4 2 2 V2 2 V4 

The parabolic fixed-point corresponds to a cusp at z = joo that is invariant 
under the parabolic elements 

J and (J i 

The volume element deriving from the hyperbolic metric reads 

dxodxidy 
dfj, = ^ , 



such that the volume of the non-compact orbifold T\H is finite. 

volfrVW) = - 0.3053218647..., 

where 

!/6^[i]-{0} 

is the Dedekind zeta function. 

We are interested in the square-integrable eigenfunctions of the Laplacian, 




that determine the quantum mechanics of a point particle moving freely in the 
orbifold T\H. We identify the solutions with Maass waveforms [18]. 

Since Maass waveforms are automorphic, and therefore periodic in xq and 
Xi, it follows that they can be expanded into a Fourier series, 

f(z) = u{y) + J2 apvKir(2*\l3\v) e 2 ^, 

/36Z[i]-{0} 

where 

[hy + hyfoy ifr = 0. 

K- lr (x) is the K-Bcsscl function whose order is connected with the eigenvalue A 
by 

A = r 2 + 1. 
If a Maass waveform vanishes in the cusp, 

lim f(z) = 0, 

it is called a Maass cusp form. Maass cusp forms are square integrable over the 
fundamental domain, (/, /) < oo, where 

(f, 9) = I fg dpi. 
Jr\n 

is the Petersson scalar product. 

According to the Roelcke-Selberg spectral resolution of the Laplacian [22], 
its spectrum contains both a discrete and a continuous part. The discrete part 
is spanned by the constant eigenfunction /o and a countable number of Maass 
cusp forms /i, /2) /3j • • • which we take to be ordered with increasing eigenvalues, 
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= A < Ai < A2 < A3 < . . .. The continuous part of the spectrum A > 1 is 
spanned by the Eisenstein series E(z, 1 + ir) which are known analytically [11]. 
The Fourier coefficients of the functions Ak (1 + w)E[z, 1 + ir) are given by 

b a = A K (1 + ir), bi = Aif(l-ir), a p = 2 ^ |-| lr , 

A, A ,ez[i] ^ 

where 



A K (s) =4ir- s T( S )( K ( S ) 

has an analytic continuation into the complex plane except for a pole at s = 1. 
Normalising the Maass cusp forms according to 

(/raj In) lj 

we can expand any square integrablc function <f> € L 2 {T\H) in terms of Maass 
waveforms, 

^ Z ) = Y / (fn^)fn(z) + ^- f (E(;s),(t>)E(z,s)ds. 

The eigenvalues and their associated Maass cusp forms are not known an- 
alytically. Thus, one has to approximate them numerically By making use of 
the Hcckc operators and the multiplicative relations among the coefficients, Steil 
[25] obtained a non-linear system of equations which allowed him to compute 
2545 consecutive eigenvalues. In [28], we extend these computations with the 
use of a variant of Hejhal's algorithm [14]. Our favourite explanation of the 
algorithm is published in [27]. 



6 Results 

The modular surface, i.e. the two-dimensional hyperbolic orbifold that is gener- 
ated by the modular group, has a reflection symmetry. This reflection symmetry 
commutes with the Laplacian. Consequently, the eigenfunctions fall into two 
symmetry classes. We call an eigenfunction to be even or odd depending on 
whether f(—x + iy) = f(x + iy) or f(—x + \y) = — fix + iy) holds. We also call 
the corresponding eigenvalue to be even or odd, respectively. 

In table 1, the first ten consecutive even and the first ten consecutive odd 
eigenvalues of the Laplacian on the modular surface are listed. 

The Picard surface, i.e. the three-dimensional hyperbolic orbifold that is 
generated by the Picard group, has two symmetries such that the eigenfunctions 
of the Laplacian fall into four symmetry classes. The eigenfunctions and their 
corresponding eigenvalues are called to be of symmetry class D, G, C, or H 
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Table 1: The first ten even and the first ten odd eigenvalues of the negative 
Laplacian on the modular surface. Listed is r. related to the eigenvalue via 



A 



odd 



13.77975135189 
17.73856338106 
19.42348147083 
21.31579594020 
22.78590849419 
24.11235272984 
25.82624371271 
26.15208544922 
27.33270808315 
28.53074769292 



9.53369526135 
12.17300832468 
14.35850951826 
16.13807317152 
16.64425920190 
18.18091783453 
19.48471385474 
20.10669468255 
21.47905754475 
22.19467397757 



Table 2: The first few eigenvalues of the negative Laplacian on the Picard 



Listed is r, 


related to the eij 


>cnvalues via A = r 2 - 


h 1. 


D 


G 


c 


H 


8.55525104 




6.62211934 




11.10856737 




10.18079978 




12.86991062 




12.11527484 


12.11527484 


14.07966049 




12.87936900 




15.34827764 




14.14833073 




15.89184204 




14.95244267 


14.95244267 


17.33640443 




16.20759420 




17.45131992 


17.45131992 


16.99496892 


16.99496892 


17.77664065 




17.86305643 


17.86305643 


19.06739052 




18.24391070 




19.22290266 




18.83298996 




19.41119126 




19.43054310 


19.43054310 


20.00754583 




20.30030720 


20.30030720 


20.70798880 


20.70798880 


20.60686743 




20.81526852 




21.37966055 


21.37966055 


21.42887079 




21.44245892 




22.12230276 




21.83248972 


21.83248972 


22.63055256 




22.58475297 


22.58475297 


22.96230105 


22.96230105 


22.85429195 




23.49617692 




23.49768305 


23.49768305 
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depending on whether 



f(x+jy) = /(iz+j) = f(-x+]y), 
f{x+iy) = /(iz+j) = -f(-x+]y), 

f(x+jy) = -ffa+i) = f(-x+jy), 
f(x + jy) = -/(ia + j) = -f(-x + jy) 

holds, respectively. 

The first few consecutive eigenvalues of each symmetry class of the Laplacian 
on the Picard surface are listed in table 2. 

Examining the eigenvalues, we find that there occur degenerate eigenvalues. 
The smallest degenerate eigenvalue of the Laplacian on the Picard surface is 
A = 12.11527484 2 + 1. The degeneracies result from the symmetries of the 
Picard surface. If desymmetrised, i.e. within a given symmetry class, there do 
not exist any degenerate eigenvalues. An explanation is given in [28]. 

Consider the level counting function 

N(r) := #{ n | n < r} 

and split it into two parts 

N(r) = N(r)+N fluc (r). 

Here N is a smooth function describing the average increase in the number of 
levels and Nfi uc describes the fluctuations around the mean such that 



1 f R 

lim - / N fluc (r)dr = 0. 



According to Wcyl's law and higher order corrections for the Picard surface 
found by Matthies [19], the average increase in the number of levels is given by 

N(r) = ) r 3 + a 2 r log r + a 3 r + a 4 

with the constants 

n — 3 



i[i| log2 + | log^ - logr(i) + § log(2 + V3) + §], 



a 4 - -j- 

If desymmetrised into the four symmetry classes, the average increase in the 
number of levels is given by 

N(r) = l^^r 3 + hr 2 + b 2 r logr + b 3 r + b A 

with the constants depending on the symmetry class as listed in table 3. 
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Table 3: The constants for the higher order corrections to Weyl's law for each 
of the four symmetry classes of the Picard surface. The constants b\ and 62 are 
known analytically [19], whereas the constants 63 and 64 have been approximated 
numerically [28]. 







b-2 


h 


b 4 


D 


1 

24 


13 
8tt 


0.8639... 


-0.288 


G 


1 

24 


3 

8tt 


0.0285... 


-0.184 


C 


1 

96 


1 

8tt 


0.0150... 


-0.062 


H 


1 

96 


1 

Sir 


0.0702... 


0.034 





2 8 4 2 s 4 



Figure 1: Level spacing distributions for the symmetry classes D (top left), 
G (top right), C (bottom left), and H (bottom right) of the Picard surface. 
The abscissae display the spacings s. The small crosses are the histograms of 
altogether 13949 nearest-neighbour spacings. The dashed curves starting at the 
origins are the integrated distributions. For comparison, the full curves show 
the Poisson distribution. 
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Unfolding the spectrum, 

Xi := N(ri), 

we are able to examine the distribution of the eigenvalues on the scale of the 
mean level spacings. Defining the sequence of nearest-neighbour level spacings 
with mean value 1 as i — > oo, 

&i ■ Xi, 

we find that the spacing distribution comes close to that of a Poisson random 
process, 

^Poisson (s) 6 , 

see figure 1, in accordance with the conjecture of arithmetic quantum chaos. 

7 Cosmology 

In the remaining sections we apply the eigenvalues and eigenfunctions of the 
Laplacian to a perturbed Robertson- Walker universe and compute the temper- 
ature fluctuations in the cosmic microwave background (CMB). 

The CMB is a relic from the primeval fireball of the early universe. It is 
the light that comes from the time when the universe was 379 000 years old. It 
was predicted by Gamow in 1948 and explained in detail by Peebles [21]. In 
1978, Penzias and Wilson won the Nobel Prize of Physics for first measuring 
the CMB at a wavelength of 7.35 cm. Within the resolution of their experiment 
they found the CMB to be completely isotropic over the whole sky. Later with 
the much better resolution of the NASA satellite mission Cosmic Background 
Explorer (COBE), Smoot et al. [23] found fluctuations in the CMB which are of 
amplitude 10 -5 relative to the mean background temperature of T = 2.725 K, 
except for the large dipole moment, see figure 2, resulting in the Nobel Prize 
for Mather and Smoot in 2006. The small fluctuations in the CMB serve as a 
fingerprint of the early universe, since the temperature fluctuations are related to 
the density fluctuations at the time of last scattering. They show how isotropic 
the universe was at early times. In the inflationary scenario the fluctuations 
originate from quantum fluctuations which are inflated to macroscopic scales. 
Due to gravitational instabilities the fluctuations grow steadily and give rise to 
the formation of stars and galaxies. 

The theoretical framework in which the CMB and its fluctuations are ex- 
plained is Einstein's general theory of relativity. Thereby a homogeneous and 
isotropic background given by a Robertson- Walker universe [12, 13, 16] is per- 
turbed. The time-evolution of the perturbations can be computed in the frame- 
work of linear perturbation theory [5] . 

An explanation for the presence of the CMB is the following, see also figure 
3: We live in an expanding universe. At early enough times the universe was so 
hot and dense that it was filled with a hot plasma consisting of ionised atoms, 
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Figure 2: Sky maps of the temperature fluctuations in the CMB as observed by 
the NASA satellite mission COBE. The sky map on the top shows the dipolc 
anisotropy after the mean background temperature of To = 2.725 K has been 
subtracted. The amplitude of the dipole anisotropy is about 3 mK. Also sub- 
tracting the dipole yields the sky map in the middle. One sees the small tem- 
perature fluctuations whose amplitude is roughly 30 )jlK. But one also sees a lot 
of foreground contamination along the equator that comes from nearby stars in 
our galaxy. After removing the foregrounds one finally gets the sky map on the 
bottom showing the temperature fluctuations in the CMB. Downloaded from 
[30]. 
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Figure 3: The expanding universe. At the time of last scattering occurcd a 
phase transition from an opaque to a transparent universe. 

unbounded electrons, and photons. Due to Thomson scattering of photons with 
electrons, the hot plasma was in thermal equilibrium and the mean free path of 
the photons was small, hence the universe was opaque. Due to its expansion the 
universe cooled down and became less dense. When the universe was around 
379 000 years old, its temperature T has dropped down to approximately 3000 K. 
At this time, called the time of last scattering, the electrons got bound to the 
nuclei forming a gas of neutral atoms, mainly hydrogen and helium, and the 
universe became transparent. Since this time the photons travel freely on their 
geodesies through the universe. At the time of last scattering the photons had an 
energy distribution according to a Planck spectrum with temperature of nearly 
3000 K. The further expansion of the universe rcdshiftcd the photons such that 
they nowadays have an energy distribution according to a Planck spectrum with 
temperature of To = 2.725 K. This is what we observe as the CMB. 

Due to the thermal equilibrium before the time of last scattering the CMB is 
nearly perfectly isotropic, but small density fluctuations lead to small tempera- 
ture fluctuations. The reason for the small temperature fluctuations comes from 
a variety of effects. The most dominant effects are the gravitational redshift that 
is larger in the directions of overdense regions, a small time delay in the tran- 
sition from opaque to transparent that slightly reduces the Hubble redshift in 
the directions of overdense regions, the intrinsic temperature fluctuations, and 
the Doppler effect due to the velocity of the plasma. 
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8 Robertson- Walker universes 

Assuming a universe whose spatial part is locally homogeneous and isotropic, 
its metric is given by the Robertson- Walker metric, 

ds 2 = dt 2 — a 2 (t)^fijdx 1 'dx\ 

where we use the Einstein summation convention. Notice that we have changed 
the notation. Instead of the quaternion z for the spatial variables, we now write 
x = xq + ixi + jx2- Jij is the metric of a homogeneous and isotropic three- 
dimensional space, and the units are rescaled such that the speed of light is 
c = 1. Introducing the conformal time dr\ = g^|y we have 

ds 2 = a 2 (rf) [drj 2 — r )ijdx l dx^ 1 \ , 

where a(rj) = d{t{rf)) is the cosmic scale factor. 

With the Robertson- Walker metric the Einstein equations simplify to the 
Friedmann equations [12, 13, 16]. One of the two Friedmann equations reads 

« + Ka 2 = ^T V + ^Aa 4 

and the other Friedmann equation is equivalent to local energy conservation, a' 
is the derivative of the cosmic scale factor with respect to the conformal time 
77. k is the curvature parameter which we choose to be negative, k = — 1. G is 
Newton's gravitational constant, Tj? is the energy-momentum tensor, and A is 
the cosmological constant. 

Assuming the energy and matter in the universe to be a perfect fluid con- 
sisting of radiation, non-relativistic matter, and a cosmological constant, the 
time-time component of the energy-momentum tensor reads 

r ° = £r (r / )+e m (ry), 

where the energy densities of radiation and matter scale like 

£r(f?) = £r(vo)(^rr-) 4 and s m (rj) = s m (j] Q )(^j^-) 3 . 

a(ri) a(rj) 

Here 770 denotes the conformal time at the present epoch. 

Specifying the initial conditions (Big Bang!) a(0) = 0, a'(0) > 0, the Fried- 
mann equation can be solved analytically [3] , 

where V{r}) = V (r\\ (72, 33) denotes the Weierstrass P-function. The so-called 
invariants gi and 53 are determined by the cosmological parameters, 

_ n A Q T 1 1 n A n r 1 n A f^ 1 

92 ~ + 12' 93 ~ ~ 16^1 216' 
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with 



_ 8nGe T (r) Q ) _ 8nGe m (r] Q ) 

where 

H(rj) 

is the Hubble parameter. 



a'(v) 
a 2 (rj) 



9 Perturbed Robertson- Walker universes 

The idealisation to an exact homogeneous and isotropic universe was essential 
to derive the spacetime of the Robertson- Walker universe. But obviously, we 
do not live in a universe that is perfectly homogeneous and isotropic. We see 
individual stars, galaxies, and in between large empty space. Knowing the 
spacetime of the Robertson- Walker universe, we can study small perturbations 
around the homogeneous and isotropic background. Since the amplitude of the 
large scale fluctuations in the universe is of relative size 10 -5 , we can use linear 
perturbation theory. In longitudinal gauge the most general scalar perturbation 
of the Robertson- Walker metric reads 

ds 2 = a 2 {ri) [(1 + 2<£)dri 2 - (1 - 2V)i ij dx i l dx j ] , 

where $ = $(77, x) and 'J = ^(rj, x) are functions of spacetime. 

Assuming that the energy and matter density in the universe can be de- 
scribed by a perfect fluid, consisting of radiation, non-relativistic matter, and a 
cosmological constant, and neglecting possible entropy perturbations, the Ein- 
stein equations reduce in first order perturbation theory [20] to 

$ = 

+ 3#(1 + c 2 )$' - c 2 A$ + (2H' + (1 + 3c 2 )(£T 2 + 1))$ = 0, 

where H = ^- and c 2 = (3 + | — ) _1 are given by the solution of the non- 
perturbed Robertson- Walker universe. In the partial differential equation for 
$ the Laplacian occurs. If the initial and the boundary conditions of 3? are 
specified, the time-evolution of the metric perturbations can be computed. 
With the separation ansatz 

^(rj,x) = J2f k (r])ip k (x)+ / dkf k (ri)il) k {x), 
k J 

where the tpk are the eigenfunctions of the negative Laplacian, and the E k are 
the corresponding eigenvalues, 

-Aip k (x) = E k ip k (x), 
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the partial differential equation for <f> simplifies to 

f' k \r,) + 3H(1 + c 2 )f k ( V ) + (c 2 s E k + 2H' + (1 + 3c 2 )(H 2 + l))/ fc (r?) = 0. 

These ordinary differential equations (one for each eigenvalue E k ) can be com- 
puted numerically in a straightforward way, and we finally obtain the metric 
of the whole universe. This gives the input to the Sachs- Wolfe formula which 
connects the metric perturbations with the temperature fluctuations, 

ST 3 f n ° d 

— (h) = 2<P(r, SLS ,x( VsLS )) - -$(0, x(0)) + 2 d V —*(»/, x(r,)), 

T 2 J VSLS drj 

where h is a unit vector in the direction of the observed photons. x(r)) is the 
geodesic along which the light travels from the surface of last scattering (SLS) 
towards us, and ?7SLS is the time of last scattering. 

If we choose the topology of the universe to be the Picard surface, we can 
use the Maass waveforms of sections 5 and 6 in the separation ansatz for the 
metric perturbations <E>. Let us further choose the initial conditions to be 

fk 0) = - = and / fc (0 = — r , 

\J E k y/E k - I 16(fi c fi r )3 

[3], which carry over to a Harrison-Zel'dovich spectrum, a is a constant in- 
dependent of k which is fitted to the amplitude of the observed temperature 
fluctuations. The quantities a k arc normal distributed random numbers. 

The following cosmological parameters are used f2 m = 0.3, 51a = 0.65, Sl c = 
1-fitot = l-nr-fim-fiAj-fffao) = 100/iokms" 1 Mpc" 1 with h = 0.65. The 
density f2 r f=a I0 -4 is determined by the current temperature To = 2.725 K. The 
point of the observer is chosen to be at a; bs = 0.2 + O.Ii + 1 . 6 j . For numerical 
reasons the infinite spectrum is cut such that only the discrete eigenvalues with 
E k = k 2 + I < 19601 and their corresponding cigenfunctions are taken into 
account. The necessary computations are carried out in [1]. The resulting sky 
map is shown in figure 4. 

Concerning the topology of the universe which manifests itself in the sup- 
pression of power in the large scale anisotropics, there exist the cosmological 
observations from COBE and WMAP. In order to quantitatively compare our 
results with these observations, we introduce the two-point correlation function, 

CV>) = (ST(n)6T(h')) coBe ^ hl . 

Figure 5 shows the correlation function corresponding to the calculated sky 
map of the Picard surface [2] in comparison with the results of the cosmological 
observations [24] and with the concordance model [24] . 

We find quite good agreement of our calculated temperature fluctuations 
with the cosmological observations, whereas the results of the concordance 
model is not in good agreement with the data for d > 7°. Especially for large 
angular separations, $ > 160°, the concordance model is not able to describe the 
observed anticorrelation in C(i9). This anticorrelation constitutes a fingerprint 
in the CMB that favours a non-trivial topology for the universe. 
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figure is taken from [1]). 
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Figure 5: Correlation functions of the calculated temperature fluctuations on 
the Picard surface (dashed line) [2] , of the WMAP observations (solid line) [24] , 
and of the concordance model (dotted line) [24] . (The figure is borrowed from 
[2])- 
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